Hierarchical cluster analysis and nonlinear mixed-effects modelling for candidate biomarker detection in preclinical models of cancer.

BACKGROUND: Preclinical models of cancer can be of translational benefit when assessing how different biomarkers are regulated in response to particular treatments. Detection of molecular biomarkers in preclinical models of cancer is difficult due inter-animal variability in responses, combined with limited accessibility of longitudinal data. METHODS: Nonlinear mixed-effects modelling (NLME) was used to analyse tumour growth data based on expected tumour growth rates observed 7 days after initial doses (DD7) of Radiotherapy (RT) and Combination of RT with DNA Damage Response Inhibitors (DDRi). Cox regression was performed to confirm an association between DD7 and survival. Hierarchical Cluster Analysis (HCA) was then used to identify candidate biomarkers impacting responses to RT and RT/DDRi and these were validated using NLME. RESULTS: Cox regression confirmed significant associations between DD7 and survival. HCA of RT treated samples, combined with NLME confirmed significant associations between DD7 and Cluster specific CD8+ Ki67 MFI, as well as DD7 and cluster specific Natural Killer cell density in RT treated mice. CONCLUSION: Application of NLME, as well as HCA of candidate biomarkers may provide additional avenues to assess the effect of RT in MC38 syngeneic tumour models. Additional studies would need to be conducted to confirm association between DD7 and biomarkers in RT/DDRi treated mice.

Preclinical Characterization of AZD9574, a Blood-Brain Barrier Penetrant Inhibitor of PARP1.

PURPOSE: We evaluated the properties and activity of AZD9574, a blood-brain barrier (BBB) penetrant selective inhibitor of PARP1, and assessed its efficacy and safety alone and in combination with temozolomide (TMZ) in preclinical models. EXPERIMENTAL DESIGN: AZD9574 was interrogated in vitro for selectivity, PARylation inhibition, PARP-DNA trapping, the ability to cross the BBB, and the potential to inhibit cancer cell proliferation. In vivo efficacy was determined using subcutaneous as well as intracranial mouse xenograft models. Mouse, rat, and monkey were used to assess AZD9574 BBB penetration and rat models were used to evaluate potential hematotoxicity for AZD9574 monotherapy and the TMZ combination. RESULTS: AZD9574 demonstrated PARP1-selectivity in fluorescence anisotropy, PARylation, and PARP-DNA trapping assays and in vivo experiments demonstrated BBB penetration. AZD9574 showed potent single agent efficacy in preclinical models with homologous recombination repair deficiency in vitro and in vivo. In an O6-methylguanine-DNA methyltransferase (MGMT)-methylated orthotopic glioma model, AZD9574 in combination with TMZ was superior in extending the survival of tumor-bearing mice compared with TMZ alone. CONCLUSIONS: The combination of three key features-PARP1 selectivity, PARP1 trapping profile, and high central nervous system penetration in a single molecule-supports the development of AZD9574 as the best-in-class PARP inhibitor for the treatment of primary and secondary brain tumors. As documented by in vitro and in vivo studies, AZD9574 shows robust anticancer efficacy as a single agent as well as in combination with TMZ. AZD9574 is currently in a phase I trial (NCT05417594). See related commentary by Lynce and Lin, p. 1217.

Mixed effects modeling of radiotherapy in combination with immune checkpoint blockade or inhibitors of the DNA damage response pathway.

Dosage optimization to maximize efficacy and minimize toxicity is a potential issue when administering radiotherapy (RT) in combination with immune checkpoint blockade (ICB) or inhibitors of the DNA Damage Response Pathway (DDRi) in the clinic. Preclinical models and mathematical modeling can help identify ideal dosage schedules to observe beneficial effects of a tri-therapy. The aim of this study is to describe a mathematical model to capture the impact of RT in combination with inhibitors of the DNA Damage Response Pathway or blockade of the immune checkpoint protein - programmed death ligand 1 (PD-L1). This model describes how RT mediated activation of antigen presenting cells can induce an increase in cytolytic T cells capable of targeting tumor cells, and how combination drugs can potentiate the immune response by inhibiting the rate of T cell exhaustion. The model was fitted using preclinical data, where MC38 tumors were treated in vivo with RT alone or in combination with anti-PD-L1 as well as with either olaparib or the ataxia telangiectasia mutated (ATM) inhibitor-AZD0156. The model successfully described the observed data and goodness-of-fit, using visual predictive checks also confirmed a successful internal model validation for each treatment modality. The results demonstrated that the anti-PD-L1 effect in combination with RT was maximal in vivo and any additional benefit of DDRi at the given dosage and schedule used was undetectable. Model fit results indicated AZD0156 to be a more potent DDRi than olaparib. Simulations of alternative doses indicated that reducing efficacy of anti-PD-L1 by 68% would potentially provide evidence for a benefit of ATM inhibition in combination with ICB and increase the relative efficacy of tri-therapy.

Radiation in Combination with Immune Checkpoint Blockade and DNA Damage Response Inhibitors in Mice: Dosage Optimization in MC38 Syngeneic Tumors via Modelling and Simulation.

Clinical trials assessing the impact of radiotherapy (RT) in combination with DNA damage response pathway inhibitors (DDRis) and/or immune checkpoint blockade are currently ongoing. However, current methods for optimizing dosage and schedule are limited. A mathematical model was developed to capture the impacts of RT in combination with DDRi and/or anti-PD-L1 [immune checkpoint inhibitor (ICI)] on tumor immune interactions in the MC38 syngeneic tumor model. The model was fitted to datasets that assessed the impact of RT in combination with the DNA protein kinase inhibitor (DNAPKi) AZD7648. The model was further fitted to datasets from studies that were used to assess both RT/ICI combinations as well as RT/ICI combinations followed by concurrent administration of the poly ADP ribose polymerase inhibitor (PARPi) olaparib. Nonlinear mixed-effects modeling was performed followed by internal validation with visual predictive checks (VPC). Simulations of alternative dosage regimens and scheduling were performed to identify optimal candidate dosage regimens of RT/DNAPKi and RT/PARPi/ICI. Model fits and VPCs confirmed a successful internal validation for both datasets and demonstrated very small differences in the median, lower, and upper percentile values of tumor diameters between RT/ICI and RT/PARPi/ICI, which indicated that the triple combination of RT/PARPi/ICI at the given dosage and schedule does not provide additional benefit compared with ICI in combination with RT. Simulation of alternative dosage regimens indicated that lowering the dosage of ICI to between 2 and 4 mg/kg could induce similar benefits to the full dosage regimen, which could be of translational benefit. SIGNIFICANCE STATEMENT: This work provides a mixed-effects model framework to quantify the effects of combination radiotherapy/DNA damage response pathway inhibitors/immune checkpoint inhibitors in preclinical tumor models and identify optimal dosage regimens, which could be of translational benefit.

An integrated modelling approach for targeted degradation: insights on optimization, data requirements and PKPD predictions from semi- or fully-mechanistic models and exact steady state solutions.

The value of an integrated mathematical modelling approach for protein degraders which combines the benefits of traditional turnover models and fully mechanistic models is presented. Firstly, we show how exact solutions of the mechanistic models of monovalent and bivalent degraders can provide insight on the role of each system parameter in driving the pharmacological response. We show how on/off binding rates and degradation rates are related to potency and maximal effect of monovalent degraders, and how such relationship can be used to suggest a compound optimization strategy. Even convoluted exact steady state solutions for bivalent degraders provide insight on the type of observations required to ensure the predictive capacity of a mechanistic approach. Specifically for PROTACs, the structure of the exact steady state solution suggests that the total remaining target at steady state, which is easily accessible experimentally, is insufficient to reconstruct the state of the whole system at equilibrium and observations on different species (such as binary/ternary complexes) are necessary. Secondly, global sensitivity analysis of fully mechanistic models for PROTACs suggests that both target and ligase baselines (actually, their ratio) are the major sources of variability in the response of non-cooperative systems, which speaks to the importance of characterizing their distribution in the target patient population. Finally, we propose a pragmatic modelling approach which incorporates the insights generated with fully mechanistic models into simpler turnover models to improve their predictive ability, hence enabling acceleration of drug discovery programs and increased probability of success in the clinic.